Accepted for publication in the Astrophysical Journal 

Preprint typeset using I^T^jX style cmulateapj v. 08/13/06 



OPTICAL SPECTROSCOPY OF BRIGHT FERMI LAT BLAZARS 

Michael S. Shaw 1 , Roger W. Romani 1 , Stephen E. Healey 1 , 
Garret Cotter 2 , Peter F. Michelson 1 , 
Anthony C. S. Readhead'* 

Accepted for publication in the Astrophysical Journal 

ABSTRACT 

We report on HET and Palomar 5 m spectroscopy of recently identified 7-ray blazars in the Fermi 
LAT Bright Source List. These data provide identifications for 10 newly discovered 7-ray flat spectrum 
radio quasars (FSRQ) and six new BL Lacs plus improved spectroscopy for six additional BL Lacs. We 
substantially improve the identification completeness of the bright LAT blazars and give new redshifts 
and z constraints, new estimates of the black hole masses and new measurements of the optical SED. 

Subject headings: BL Lacertae objects: general — galaxies: active — quasars: general — surveys 



1. INTRODUCTION 

The Fermi LAT pair production telescope has been 
surveying the 20 MeV-300 GeV 7-ray sky since 2008 Au- 
gust 11. Among the many sources being detected by this 
mission, blazars dominate the extragalactic sky. The first 
published set of these objects, the 'LAT Bright AGN 
Sample' (LBAS) based on 3 months of sky survey expo- 
sure, included 132 bright (> 10cr) detections at \b\ > 10°. 
The bulk (117) of these sources have been associated 
with flat spectrum radio counterparts, which in turn end 
up being mostly well-known flat spectrum radio quasars 
(FSRQ) and BL Lac objects. A few additional sources 
are pulsars and a few remain unidentified. The associ- 
ations are based on the distribution of radio and X-ray 
properties of the radio counterparts (CRATES; Healey 
et al. 2007) and the analysis allows a quantitative as- 
sessment of the probability of association. 

Even before the start of the Fermi mission, it was 
recognized that a large sample of these blazars would 
be needed for source identification; in the Candidate 
Gamma-Ray Blazar Survey (CGRaBS; Healey et al. 
2008, and references therein) we pursued optical spec- 
troscopy of the sources selected to be most similar to 
the previously known EGRET blazars. This set of 1625 
blazars included many sources without previous spectro- 
scopic identification. Here we report on previously un- 
published optical spectra from our CRATES/CGRaBS 
survey and observations for new sources from the LBAS 
list. The data allow classification of the blazar type and 
redshift solutions for the strong-lined objects. Significant 
rcdshift constraints are also obtained for the weak-lined 
BL Lacs. These measurements contribute ten FSRQ 
redshifts, one new spectroscopic redshift for a BL Lac, 
and 11 additional BL Lac redshift constraints. Includ- 
ing these new data, 91% of all LBAS sources have been 
spectroscopically typed, 84% have redshifts, and 70% of 
the BL Lac have spectroscopic redshifts. We also use 
our spectra to examine the black hole masses and optical 
continua of these blazar sources. 

1 Department of Physics/KIPAC, Stanford University, Stanford, 
CA 94305 

2 Department of Astrophysics, University of Oxford, Oxford 
OX1 3RH, UK 

3 Department of Astronomy, California Institute of Technology, 
Pasadena, CA 91125 



In this paper, we assume an approximate concordance 
cosmology — £l m — 0.3, f^A = 0.7, and H = 70km s _1 
Mpc" 1 . 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. HET 

The bulk of the spectroscopic observations were per- 
formed with the 9.2 m Hobby-Eberly Telescope (HET) 
at McDonald Observatory; many o f these observation s 
were part of the CGRaBS survey (|Healev et al.l 12008). 
but additional data were taken on newly discovered 
LBAS sources. The HET observes in the declination 
range —11° < i5 < + 73°, and is fully queue scheduled 
l|Shetrone et al.l l2007t ), allowing us to receive data re- 
motely year round and to spread the cost of inclement 
weather and unfavorable conditions among the observing 
programs. We use the Marcar io Low-Resolution Spectro- 
graph, LRS (|Hill et al.lll998l ). with grism Gl (300 lines 
mm -1 ), a 2" slit, and a Schott GG385 long-pass filter 
for a resolution of R w 500 between 4150 A and 10500 
A. Typical exposures are 2 x 600 s for FSRQ objects and 
2 x 900 s for BL Lac objects, with the slit placed along 
the parallactic angle. 

2.2. DBSP 

We also used the double spectrograph (DBSP) on the 
5 m Hale Telescope at Palomar. The observations with 
a 1" slit at the parallactic angle used a 600 line mm -1 
grating on the blue side, covering AA3100 — 5200 at a 
resolution of - 2.8 A (R ~ 1450). The red camera, 
with a 316 line mm -1 grating covers AA5200 — 9500, at 
a resolution ~ 5.2 A (R ~ 1350). P200/DBSP is our 
primary tool for targets outside of the HET declination 
range (i.e. -35° < 5 < -11° and +73° < S) and provides 
extra spectral coverage and resolution in the blue for 
bright BL Lacs. Typical exposures were 2 x 600 s for 
FSRQs and up to 3 x 900 s for faint BL Lacs. 

2.3. Analysis Pipeline 

Data redu ction was performed with the IRAF package 
(|Todvlll986r ) using standard techniques. Wavelength cal- 
ibration was performed with a neon-argon lamp at the 
HET and on the red side at DBSP, and iron-argon on 
the DBSP blue side. For these relatively faint objects, 



2 



M. S. Shaw et al 



3 - 



1 





2 

1.5 
1 

0.5 




1 - 



O 

X 

X 

p 
U- 

%- 
o 

0) 

Oh 



0.5 


0.3 

0.2 

0.1 


0.3 



0.2 - 



0.1 






4000 



J0017-0512 




1.5 



0.5 




6000 8000 4000 6000 8000 

Wavelength (A) Wavelength (A) 

Fig. 1. — New spectra of FSRQ objects in the Fermi LBAS sample. The emission line properties are listed in Table 1. 



we employ an optimal extraction algorithm ( Valde3 fl992f l 
to maximize the S/N. 

We perform s pectrophoto metric calibration using stan- 
dard stars from I Oka (|1990f K In most cases the standard 
exposures from the data night were suitable, but at the 
HET, standards from subsequent nights were sometimes 
used. Due to differential slit losses and variable condi- 
tions between object and standard star exposures, we 
estimate that the accuracy of our ab solute spectropho- 
tometry is ~ 30% (jHealev et alJl2008t ). although the rel- 
ative spectrophotometry is considerably better. Spectra 
are corrected for telluric absorptions and visually cleaned 
of cosmic rays. Multiple exposures on a single target are 



combined into a single spectrum, weighting by S/N. 

We visually identify and measure emission line equiva- 
lent widths and FWHMs and derive redshifts from these 
line measurements. They are listed in Table 1. Line and 
redshift measurement techniques are discussed in §3.1. 



3. RESULTS 

We present spectra for ten FSRQ objects, six new 
BL Lac identifications and spectra, and high S/N spec- 
troscopy on six additional objects previously known to 
be BL Lacs. We present spectroscopic redshifts for all 
ten FSRQ, and for one BL Lac. We derive virial BH 
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TABLE 1 

Line Properties for New LBAS FSRQ 
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masses for the FSRQ, and for BL Lacs that lack spectro- 
scopic redshifts, we extract significant constraints on the 
redshift. 

3.1. FSRQ Spectra 

We present ten previously unpublished FSRQ spectra 
with redshifts ranging from 0.2 27 to 1.805. Thes e red- 
shifts were previously listed in lAbdo et al.1 j2009:). All 
redshifts are confirmed by multiple emission lines and dc- 
rive d by cross-correlation analysis using the rvsao pack- 
age (|Kurtz fe M ink 1998). Table 1 gives the approximate 
continuum fluxes at 10 147 Hz (A = 5981 A) and spectral 
indices a (F v cx v~ a ). The fluxes, as above, have un- 
certainties as large as 30% due to unknown slit losses. 
The spectral index uncertainties are estimated by fits to 
independent subsets of the spectral range. 

The properties of strong emission lines are listed in 
Table 1. We are particularly interested in the broad 
emission lines H/3, Mg II, and C IV, as their virialized 
broad components are used to make estimates of the 
black hole mass. For these lines, we take care to measure 
the Gaussian full width at half maximum (FWHM) of 
the broad component, avoiding contamination by other 
spectral features. 

For H/3, it is important to include narrow components 
in the lin e fit. We fit the cont i nuum to a power law, and, 
following iMcLure fc Dunlopl (|2004l ), we simultaneously 
fit broad and narrow H/3, narrow [O IIIJ4959 A, and nar- 
row [O III] 5007 A. We require rest wavelengths of narrow 
H/3 and the [O III] lines to match laboratory values; the 
broad H/3 wavelength was free. All lines are modeled 
with Gaussian profiles. The continuum is measured at 
5100 A. 

For Mg II, we first sub tract off our power l aw fit. Then, 
we subtract a template (|Tsuzuki et al|[2006l ) of the broad 
Fe II complexes near 2800 A, to minimize contamina- 
tion with Mg II. We then fit the Mg II line itself with 
a broad+narrow Gaussian. For this line the continuum 
measurement is , as usual, made at 3000 A, to avoid Fe 
contamination (|McLure fc Dunlodl200l . 

C IV lines sometimes suffer strong associated absorp- 
tion, but such corrections were modest for the few cases 
presented here. Also, for consistency with the SDSS mea- 
surements (see below) and to minimize sensitivity to as- 
sociated absor bers we report he re the FHWM of the ac- 
tual C IV line (IShen et al.1 12001 . 



3.2. Black Hole Mass Estimates 

We can use the emission line measurements reported 
above to estimate black hole masses for the FSRQ ob- 
jects in our s ample using the empir ical 'viriaP scaling 
relationships (|McLure fc Jarvisl 120021 ). Mass scaling re- 
lations have been applied to large samples of quasars in 
the past a nd are calibrated by reverb eration mapping 
techniques (jVestergaard fc Osmerll2009t ). 




0.5 1 1.5 2 2.5 

redshift z 

Fig. 2. — BH mass estimates as a function of redshift for LBAS 
sources. Filled points are measured from FSRQ spectra presented 
in this paper, open points are for spectra of LBAS FSRQ available 
from the SDSS archive (Squares: Hf3 A4861, Triangles: Mg II 
A2800, Circles: C IV A1550). When two species are used for a given 
source, dashed lines connect the mass estimates. For comparison, 
the mea n and rms range of mass estimates from the full SDSS 
catalog l lShen 'eFaTl 2008) are shown. 

Continuum flux measurements were made for each line 
at the appropriate rest wavelength (5100 A for H/3, 3000 
A for Mg II, 1350 A for C IV), and converted to contin- 
uum luminosity \L\, using our assumed cosmology. For 
two objects in our sample, J1012+2439 and J2327+0940, 
the C IV line is clearly present, but the spectra do not 
extend blue-ward enough to cover 1350 A; for these, we 
extrapolate a power law fit to the continuum at longer 
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wavelengths. We follow IShen et all (|2008l ) in the calcu- 
lation of the black hole masses: 

f M \ , , / \L\ \ , / F WHM \ 

(i) 

where a = 0.505, b = 0.62 for Mg II; a = 0.672, b = 0.61 
for H/3; and a = 0.66, b = 0.53 for C IV. 

Masses are plotted in Figure [2] Values measured from 
spectra presented in this paper are indicated by filled 
points (species used indicated by point style). Some 
of the LBAS FSRQ h ave spectra in the SDSS archive 
(Schneider et all 120071 ) . We recovered these and mea- 
sured line profiles and continuum fluxes exactly as for 
our new objects to obtain similar BH mass estimates 
(open points). For about half of these sources, SPSS pub- 
lished independently estimated values for M (|Shen et al.l 
2008). We find adequate agreement, with rms differences 
in log(M) of 0.04 (H/3) and 0.14 (Mg II, C IV), compa- 
rable to the inter-species dispersion. 

In general, the virial mass estimators are of limited ac- 
curacy for individual objects. As one illustration, when 
estimates are available for a given object from more than 
one species, we connect the points with a dashed line. 
The pipeline processing of the entire Sloan DR5 sample 
of optically selected quasars provides a large number of 
such mass estimates. In Figure [2] we show the mean and 
rms log(M) values as a function of redshift. For both 
the 7-ray selected sample presented here and the opti- 
cal SDSS sample, the mean mass increases with redshift, 
likely a simple flux bias. However, we note that the 7-ray 
sources trend about 3x lower in mass than the optically 
selected objects, albeit with large scatter. There are sev- 
eral possible origins of this trend. One possibility is that 
this is a result of the fainter optical limit in the follow-up 
to the LBAS sample. Another selection effect may be im- 
posed by the 7-ray selection's dependence on the jet flux, 
which, as a strong function of the viewing angle to the 
relativistic jet, is more weakly dependent on the underly- 
ing source luminosity and presumably black hole mass. A 
final possibility stems from orientation: virial estimates 
(Equation 1) assume isotropic velocities in the broad line 
region. However, blazars are known to have the jet axis, 
and hence the black hole spin and inner accretion disk 
axis, close to the Earth line-of- sight. If the bro ad line 
region has a toroidal structure ([Fine et al.l 12008). then 
blazars will have lines with FWHM decreased by (sini) 
and lower virial mass estimates than QSOs viewed at 
larger inclination angles. Larger data sets from future 
LAT samples should allow us to probe the reality and 
origin of such effects. 

3.3. BL Lac Properties 

We have obtained high sensitivity (S/N ~ 100) spec- 
tra of twel ve LBAS BL Lacs. Following iMarcha et all 
(1996) and lHealev et al.1 (|2008l ). we adopt here the prag- 
matic 'optical spectroscopic' definition of a BL Lac as 
a blazar lacking emission lines with observed equivalent 
width greater than 5 A and a limit on any possible 4000 
A spectral break of < 40%. For some sources in CGRaBS 
we have seen that, as the continuum level fluctuates, the 
source passes from spectroscopic BL Lac to spectroscopic 
FSRQ. We designate a source as a BL Lac if during any 
of our observations with adequate spectral range it has 



satisfied the spectroscopic definition. In this connection, 
one of the high-mass objects appearing in Figure[2]should 
be mentioned. This is J1058+0133, designated an FSRQ 
in LBAS based on a clear detection of a broad Mg II 
line at z=0.888. However, this line has an EW of only 
3.5 A; accordingly this blazar should be formally des- 
ignated a BL Lac at the epoch of the SDSS spectrum. 
The very strong non-thermal continuum, in fact, drives 
up the apparent luminosity and the BH mass estimate. 
This emphasizes that non-thermal continuum can arti- 
ficially increase the apparent mass; the virial estimates 
are not calibrated for BL Lac sources. 

The spectral range requirement for BL Lac designa- 
tion is that the spectrum covers [A mm to Aniax] with 
sufficient S/N to detect at least one of the standard 
AGN broad emission lines — Lya, C IV, C III, Mg II, 
H/3, or Ha. In practice this means adequate S/N over 
Amax/A mm > 1.75; this is easily satisfied for the BL 
Lacs in our LBAS sample. 

To describe the BL Lac continua, we fit the spectra 
to a power law, using the IRAF nfitld routine. Most 
are well fit with rising (blue) power-laws, suggesting 
a synchrotron vF v peak well above the optical band. 
J0144+2705 is, however, very red. This is confirmed 
by USNO B1/2MASS fluxes which rise to the J-band 
and flatten above 1.5/xm. The Galactic extinction in the 
direction of this source is Ab = 0.3 mag, so local absorp- 
tion is unlikely to redden the source. It seems plausible 
that this BL Lac has a low synchrotron peak frequency 
z'max ~ 10 14 Hz. This may be tested by broader study 
of the SED. 

The spectral indices (F v cx v~ a ) and fluxes of the BL 
Lac are reported in Table [2] For the continuum nor- 
malization we assume the 30% absolute spectrophoto- 
metric uncertainty noted above; for the spectral indices 
we again estimate the dominant systematic errors from 
fits to non-overlapping sections of the continuum. The 
continuum-normalized spectra are shown in Figures 3 & 
4. These plots also label notable interstellar features e.g. 
Na A5892, an d the diffuse interstellar bands, or DIBs 
(|Herbiglll995h . 

3.4. BL Lac Redshift Limits 

The strong dominance of the non-thermal nuclear con- 
tinuum in BL Lacs guarantees that the equivalent widths 
of any broad nuclear lines or absorption features from the 
host galaxy are very small. While recent work with 8 
m-class telescopes has provided some redshift solutions, 
spectroscopic redshifts for many objects remain elusive 
(Sbarufatti et al 20051) ■ Nevertheless, quantitative lim- 
its on host spectral features can be useful in extracting 
redshift constraints. For example, based on the lack of 
a Lya forest-induced break in the HET spectra, none of 
our BL Lacs can have z > 2.41. This rather weak bound 
can be improved with UV coverage beyond the HET LRS 
4150 A limit. 

Also, there is good evidence, especially from HST 
imaging, that BL Lacs reside in relati vely uniform gi- 
ant e lliptical hosts, with M ~ —22.9 ([Sbarufatti et all 
2005) . Analysis to limit the contribution of such a host 
to the observed spectru m can thus provide a lo wer limit 
on the BL Lac redshift ([Sbarufatti et alll2006h . We de- 
velop here a similar spectral bound on redshift, using the 
Ca H/K and G-band absorption features, the strongest 
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Fig. 3. — New spectra of LBAS BL Lac objects. Top panels are flux-calibrated spectra with flux in 10 erg cm s 
panels show spectra normalized by dividing out the continua. Telluric and interstellar absorption features are labeled. 
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narrow features in the host spectrum. 

We first determine the expected equivalent width of a 
given absorption line by measuring an elliptical template 
spectrum and then computing the net equivalent width 
in the ob served host+nucleus sp ectrum at each redshift 
following ISbarufatti et al (|2006h : 



EW, 



_ (1 + z)EW h 
expected- 1 + F n /F h 



(2) 



where EW-^ is the equivalent width of the line in the host 
galaxy spectrum, F n is the nuclear flux at the wavelength 
of that absorption line, and is the flux of the giant 
elliptical host at that wavelength. F^ is corrected for the 
slit/extraction aperture losses, assuming a deVaucoleurs 



profile with redshift-dependent angular size set by a con- 
stant radius ro = lOkpc. Fjx is the observed flux minus 
that expected from the host in our extraction aperture. 
This predicted feature strength is computed for the ma- 
jor absorption lines at each redshift. 

We next compute the 3ct limit on narrow absorption 
features in the observed spectrum as a function of wave- 
length from the local rms of negative spectral fluctua- 
tions, after high-pass filtering to remove broad contin- 
uum features and after exclusion of zero-redshift Galac- 
tic and telluric features. We find that the line sensitivity 
varies considerably across the spectrum (Sbarufatti et al. 
2006 assume constant sensitivity for their VLT spectra). 
For the HET data, limited blue sensitivity weakens the 
constraints below ~ 5000 A, while strong fringing (uncor- 
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Fig. 4. — As for Figure 3. J1054+2210 has a redshift solution z 
z =1.04. 

rectable due the varying HET pupil) weakens the bounds 
beyond ~ 8000 A. The sensitivity also varies due to the 
changing resolution; this is especially important for the 
DBSP data, although no DBSP spectra are used in the 
measurements of the BL Lacs presented here. 

Comparison of the expected and limiting EW curves 
gives a lower limit on the BL Lac redshift (for the as- 
sumed standard host magnitude) whenever one absorp- 
tion feature is excluded at the 3a level. Generally the 
Ca features provide the strongest constraint. However, 
since our HET spectra are limited to A > 4150 A, the G- 
band at A = 4304 A is used to exclude z < 0.05. Figure 
[5] shows the exclusion curve for J1253+5301, where we 
find z > 0.77. The results for the BL Lac data are pre- 
sented in Tabled With our HET data, we can exclude 



6000 8000 
Wavelength(A) 
=0.634. J2325+3957 has an intervening absorption line system at 

redshifts less than z ~ 0.5 — 0.8. For a few particularly 
bright BL Lac, however, the redshift constraints are sub- 
stantially weaker as the expected equivalent widths are 
smaller. Higher resolution, higher S/N spectra can, of 
course, improve these bounds. 

A few objects deserve particular comment. For 
J1054+2210, after completing the redshift limit analysis, 
weak H/K and G-band features plus O II A3727 emis- 
sion were found at redshift z = 0.634. This detection 
is confirmed by cross-correlation analysis. The measured 
redshift is just above the redshift bound of z > 0.60 given 
by the exclusion analysis. For J2325+3957, the EW limit 
analysis gave a redshift bound of z > 0.77. We have 
detected an intergalactic absorption system in this spec- 
trum, with Mg II doublet and Fe II absorption features 
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TABLE 2 

BL Lac Continuum properties for LBAS sources 



Fermi Name 
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RA 


dec *l,10".7 
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14 27 00.3 
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0.01±0.01 


> 0.03 


OFGL J1543.14 
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15 42 56.9 


61 29 55.3 


61 


0.74±0.02 


> 0.63 


OFGL J1719.34 


1746 


J17194-1745 


17 19 13 


17 45 06.3 


25 


0.67±0.03 


> 0.58 


OFGL J2325.3+3959 


J23254-3957 


23 25 17.8 


39 57 36.5 


17 


0.76±0.02 


> 1.05 



Note. — t Measured redshift z = 0.634. 
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Fig. 5. — Expected equivalent widths of the H (3934 A, triangles) 
and G (4304 A, squares) features calculated at different redshifts. 
The 3<t limit on line equivalent widths across the spectrum is plot- 
ted with open circles. As long as either H or G is above the circles, 
the corresponding redshift is excluded. Here, we excluded z < 0.77. 

corresponding to z = 1.04. This stronger bound, listed 
in Table [21 is consistent with the EW limit analysis. 

4. CONCLUSIONS 

The optical spectroscopy presented here provides sig- 
nificant new information on some of the less well-known 
members of the set of the LBAS 7-ray blazars. 

For all of the new FSRQ, precision redshifts are mea- 
sured from the spectra. For one of the new BL Lacs, 
we determine the redshift spectroscopically; for the oth- 
ers we extract limits on the redshift from the spectrum. 
Continuum fluxes and spectral indices are provided that 
will be useful for SED plots. These measurements bring 



the LBAS to 84% redshift completion. For most of the 
unsolved BL Lacs we extract significant lower limits on 
the redshift. These will be helpful in modeling the 7-ray 
BL Lac population and evolution and in focusing redshift 
searches to a narrower spectral range. 

For the strong-line objects we present estimates of 
the virial black hole masses. Although such mass es- 
timates have large uncertainties, our objects appear to 
have somewhat lower mass than typical for the SDSS 
quasar sample. Like that sample, the mean mass in- 
creases with redshift. In the one case where a BL Lac- like 
object has a mass measurement, the excess continuum 
artificially inflates the apparent black hole mass. 

The spectroscopic analysis presented here for the small 
LBAS blazar sample will be much more powerful when 
applied to the full sample of LAT blazars detected during 
the first year sky survey, which should include 500-1000 
objects. We continue to collect spectroscopic data on 
the LAT blazars as they are detected and plan to use 
the optical analysis of this larger data set to probe the 
evolution of the 7-ray emitting objects. 



The Hobby*Eberly Telescope (HET) is a joint project 
of the University of Texas at Austin, the Pennsyl- 
vania State University, Stanford University, Ludwig- 
Maximilians-Universitaet Muenchen, and Georg-August- 
Universitaet Goettingen. The HET is named in honor of 
its principal benefactors, William P. Hobby and Robert 
E. Eberly. The Marcario Low Resolution Spectrograph is 
named for Mike Marcario of High Lonesome Optics, who 
fabricated several optics for the instrument but died be- 
fore its completion. The LRS is a joint project of the 
Hobby*Eberly Telescope partnership and the Instituto 
de Astronoma de la Universidad Nacional Autonoma de 
Mexico. 

MSS and SEH are supported by NASA under grants 
NNX08AW30G and NAS5-00147. ACSR is supported by 
the NSF under grant AST-0808050. 
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